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ABSTRACT 

We investigate the composition of interstellar grains along the line of sight toward C Ophiuchi, a well- 
studied environment near the diffuse-dense cloud transition. A spectral decomposition analysis of the 
solid-state absorbers is performed using archival spectroscopic observations from the Spitzer Space 
Telescope and Infrared Space Observatory. We find strong evidence for the presence of sub-micron- 
sized amorphous silicate grains, principally comprised of olivine-like composition, with no convincing 
evidence of H 2 O ice mantles. However, tentative evidence for thick H 2 O ice mantles on large (a Ri 
2.8 fim) grains is presented. Solid-state abundances of elemental Mg, Si, Fe, and O are inferred from 
our analysis and compared to standard reference abundances. We find that nearly all of elemental Mg 
and Si along the line of sight are present in amorphous silicate grains, while a substantial fraction of 
elemental Fe resides in compounds other than silicates. Moreover, we find that the total abundance of 
elemental O is largely inconsistent with the adopted reference abundances, indicating that as much as 
^156 ppm of interstellar O is missing along the line of sight. After taking into account additional limits 
on the abundance of elemental O in other 0-bearing solids, we conclude that any missing reservoir of 
elemental O must reside on large grains that are nearly opaque to infrared radiation. 

Keywords: dust, extinction — local interstellar matter — infrared: ISM — ISM: abundances — ISM: 
molecules — stars: individual {( Oph, HD 149757) 


1. INTRODUCTION 

The interchange of gas and dust plays an essential 
role in the evolution of the interstellar medium (ISM). 
The surfaces of the dust grains act as both reposito¬ 
ries for atoms and molecules adsorbed from the gas- 
phase and active sites for molecule formation over a wide 
range of interstellar environments. The apparent under¬ 
abundances or “depletions” of many heavy elements de¬ 
tected in the gas-phase are interpreted as evidence for 
their inclusion on or in interstellar grains. Forty years 
ago, iGreenbei^ (II974D concluded that for the cosmi- 
cally abundant elements of carbon, nitrogen, and oxy¬ 
gen, significantly larger amounts were missing from the 
gas-phase than could be accounted for by accretion onto 
interstellar dust grains. In view of this critical finding, 
iGreenbeT^ (|I974fl hypothesized that “there is more hid¬ 
den mass than there is mass in the dust”, with such mass 
possibly existing in the form of intermediate-sized “snow¬ 
balls” or more complicated interstellar molecules that 
had yet to be detected in the solid-state. 

With major advancements in the sensitivity and 
resolution of ultraviolet spectroscopy, as well as de¬ 
creases in solar and stellar abundance estimates, ap¬ 
parent inconsistencies between the expected (approxi¬ 
mately solar) and observed (gas and dust) interstellar 
abundances have la.rgelv subsided oyer the years (e.g. , 

Gardelli et al.lll994j^Sofia et al.llI99^lSnow fc Wittl[l996t 

Sofia fc Mevedl200lll . However, the depletion of elemen¬ 

tal O remains problematic in the more heavily depleted 
regions of the ISM. In a recent comprehensive study, 
l.lenkin^ (|2009[ 1 presents a unified representation of de¬ 
pletions, in which previous results from the literature 


are re-analyzed using a consistent system of oscillator 
strengths for the various transitions used to derive gas- 
phase column den sities. Independ ent of the adopted ref¬ 
erence abundance, I.IenkinsI (j2009f l finds that elemental O 
(as observed in its dominant neutral atomic form, O I, in 
the diffuse ISM) is depleted from the gas-phase at a rate 
that far exceeds the rate at which it can be sequestered 
into silicate and metallic oxide particles. 

As later explored in detail bv iWhitt^ 1)201(1( 1. this re¬ 
sult implies a crisis in our understanding of interstel¬ 
lar dust. At the high end of the density range covered 
by the depletion measurements (i.e., near the transition 
between diffuse neutral and dense molecular phases of 
the ISM), as much as 25%-30% of the total elemental O 
budget is unaccounted for in any known observed form. 
Other than silicates and oxides, possible 0-bearing reser¬ 
voirs in com mon grain materials in clude organic refrac¬ 
tory matter l)Li fc Greenberg|ll997fl and interstellar ices 
of H 2 O, CO, and CO 2 . However, ice mantles h ave only 
been detected in the denser regions of the ISM (jWhittetl 
1201(1 and references therein), and a failure to detect the 
^6 /im C=0 carbonyl feature in diffuse ISM dust ap¬ 
pears to place strong co nstraints on the abundance of O- 
bearing organic m atter ((Pendleton fc Allamandoi8]l2002l : 
iChiar et al.ll2013ll . 

A further possibility is that 0-bearing gaseous 
molecules, such as CO and O 2 , might be important con¬ 
tributors. However, in the case of CO, the observed 
abundances toward diffuse molecular clouds are well- 
constrained, and clearly too small to account for a sig- 
nificant fraction of the total elemental O abundance 
(ILiszt et al.l 1200911 . Abundances of O 2 are less well- 
constrained, but are unlikely to exceed those of CO. Fi- 
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Figure 1. Combined archival 750-SWS and ^p^t-zer-IRS spectra 
(black lines) of the Oph line of sight and 3.4 fim broadband pho¬ 
tometry (magenta circle) from WISE^ shown with their respective 
uncertainties (gray lines). For comparison, the expected photo- 
spheric continuum emission by an 09. 5 V star (blue line), cor re- 
sponding to a T = 34300 K blackbody HHowarth &: Smithll^OOHh is 
also indicated. Note that in the Rayleigh-Jeans limit, a blackbody 
spectrum has equal to a constant. 

nally, the possible existence of large, micron-sized inter¬ 
stellar grains have also been proposed as viab le 0-bearing 
repositories (iGreenberd 11974 iJenkinsI I2009D , in view of 
the known difficulty in detecting such materials. 

In this paper, we re-examine the line of sight toward 
the bright 09.5 V star ^ Ophiuchi Oph, HD 149757; 
iHerbid [l968f) . The low-velocity {vq = —15 km s“^) 
component toward this well-studied star is regarded as 
a prototypical “cool diffuse” cloud with high elemen¬ 
tal depletions and partial H I —>■ H 2 conversion, and is 
representative o f an environment near the diffuse-dense 
ISM transition (Savage &: Sembaci] 119961 : IJenkinsI[20091 : 
iLiszt et ahl l2009|l . We utilize arch ival infrared spectra 
from the Spitzer Space Telescope (IWerner e t al.112004f) 
and t he Infrared Space Observatory f/5'0: [Kessler et al.1 
[TM^ to quantify the solid-state absorption by refrac¬ 
tory and molecular compounds along the line of sight, 
and combine our results with previously published data 
from the literature to assess the overall elemental budget 
toward C Oph. 

In general, '^20-29% of the elemental O abundance 
in the ISM is expected to be present in amorphous sil¬ 
icate grains. The chemical and physical properties of 
amorphous silicates can be determined by detailed mod¬ 
eling of their observed 10 and 20 /rm resonance fea¬ 
tures. Revealed through infrared spectroscopy, it is well- 
known that the profile of the 10 /iin sili cate emission 
featu re detected toward the Trapezium (|Forrest et al.1 
[19^ is similar to the silicate abso rption feature ob¬ 
served toward dense in olecular clouds (|Gillett et al.lll97^ 
iRoche fc Aitkeniri984ll . whereas silicates in diffuse clouds 
possess a narrower absorption profile (|Roche fc AitkenI 


11984 iBowev et aH 1199811 . However, the high signal-to- 
noise ratio of the Spitzer observations toward C, Oph al¬ 
lows us, for the first time, to characterize the composition 
of the silicates in a diffuse interstellar cloud near the H I 
—H 2 transition. 

In Section [21 an overview of the archival ISO and 
Spitzer observations is presented. In Section O spectral 
decomposition analyses of the infrared spectra toward 
C, Oph are employed to investigate the presence of solid- 
state absorption along the line of sight. In Section 01 we 
use our results to quantify the elemental abundances of 
Mg, Si, O, and Fe carried by amorphous silicate grains, 
and also place strict limits on the abundance of elemen¬ 
tal O in other 0-bearing solids. Finally, implications of 
our results with regard to the “missing oxygen” problem 
are discussed in Section [SI 

2. ARCHIVAL OBSERVATIONS 

Observations toward C, Oph are assembled using 
archival spectra obtained by the Spitzer Spaee Telescope 
Infrared Spectrograph (IRS; iHouck et al.l 120041 and the 
Short Wavelength Spectrometer fSWS: Ide Graauw et HI 
1199611 on board ISO. A s part of the Spitzer Atlas of Stel¬ 
lar Spectra (PID 485: lArdila et al.11201011 . the Spitzer- 
IRS observations {Spitzer AOR key 27582976) were car¬ 
ried out with the low-resolution {R = 60-120) mod¬ 
ules, Short-Low (SL; 5.2-14.0 ^m) and Long-Low (LL; 
14.0-36.1 /im), on 2008 September 10. The observa¬ 
tions utilized the standard two-position nod configura¬ 
tion, yielding a total on-source exposure time of 12 s for 
each module. The extracted 5-36 pLui Spitzer-IRS spec¬ 
trum was retrieved from the Cornell Atlas of Spi tzer IRS 
Sources (CASSIS v5.2: iLebouteiller et al.llMo . which 
is based on the Advanced Optimal Extraction of the 
Spitz er Science Center SI8.18 .0 pipeline basic calibrated 
data ([Lebouteiller et afll2010ll . Post-extraction process¬ 
ing was performed in order to scale the flux density of 
the SL module to that of the LL module at 14 /im. The 
uncertainties in the flux density are estimated to be half 
the difference between the two independent spectra from 
each nod position. 

Covering the nominal spectral range of 2.4-45 /iin, the 
ISO-SWS full grating scans were performed on 1997 Au¬ 
gust 5 {ISO TDT 62803102) with the SWSOl Astronom¬ 
ical Observing Template at speed 1. The observations 
result in a mean spectral resolving power oi R 200- 
300 and a total exposure time of 1140 s. The extracted 
ISO-SWS spectrum was r etrieved from the SWS Atlas 
of fully processed spectra (jSloan et al.l 120031) . Utilizing 
broadband photometr y from the Wide-fie ld Infrared Sur¬ 
vey Explorer {WISE-. I Wright et al.ll2010l) All-Sky Data 
Release, the ISO-SWS spectrum was scaled by a factor 
of 1.01 to match the flux density derived from the WISE 
3.4 /rm passband. 

The combined 2-36 /im spectra are shown in Figure 
[TJ However, due to the relatively poor signal-to-noise ra¬ 
tio in the ISO-SWS spectrum at wavelengths longward 
of 4 /im, only the 2.4-4 /im spectral range is retained 
to investigate the presence of H 2 O ice along the line of 
sight. For comparison, the expected photospheric con¬ 
tinuum emission by an 09. 5 V star, corresponding to 
a T = 34300 K blackbody (jHowarth fc Smith! 12001!) . is 
also indicated. Note that in the Rayleigh-Jeans limit, a 
blackbody spectrum has }?Fy equal to a constant. Dif- 
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Figure 2. Spitzer-IKS spectrum and best-fit spectral decomposition models for the Oph sightline. Spectral decomposition models are 
shown for a second- (left panels) and fourth-order (right panels) polynomial flux continuum; the third-order spectral decomposition model 
results are identical to those of the fourth-order spectral decomposition model. Top panels: observed spectrum (black line), model flux 
continuum (gray line), and best-fit model spectrum (blue line). Bottom panels: observed optical depth spectrum (black line) and best-fit 
model optical depth (blue line). The total absorption is modeled as a linear combination of amorphous silicates with Mg 2 Si 04 (green line), 
MgFeSi 04 (violet line), and MgSiOs (red line) stoichiometries. 


ferences between the observed spectrum and that of a 
pure stellar photosphere are attributed to extinction by 
intervening dust, consisting of amorphous silicates and 
possibly other dust species lacking spectral features in 
the mid-infrared wavelength range (e.g., amorphous car¬ 
bon or metallic iron). 

3. SPECTRAL DECOMPOSITION ANALYSIS 

In order to characterize the interstellar absorbers along 
the line of sight toward C Oph, we approximate the atten¬ 
uation by dust as an intermediate screen against a bright 
infrared continuum. Historically, infrared absorption fea¬ 
tures have been studied by modeling optical depth spec¬ 
tra with laboratory-derived spectra of candidate mate¬ 
rials. However, in this classical approach, optical depth 
spectra are often inferred through the subjective pro¬ 
cess of local-continuum fitting. In consequence, this pro¬ 
cess may artificially set the underlying optical depth to 
be zero over a pre-defined spectral region, which can 
strongly affect the subsequent modeling of the absorp¬ 
tion profiles. In light of this concern, we choose to carry 
out a different method in extracting information on the 
absorbers. 


In our approach, which is based on simulating the 
emergent spectral energy distribution (SED), the flux 
continuum and extinction opacity are simultaneously de¬ 
termined in the logarithmic domain by 

(I) 

where and E™"* are the observed and continuum 
flux densities, respectively, T\^i is the normalized optical 
depth of an individual absorbing component, and ai is 
its corresponding non-negative scaling parameter. The 
flux continuum is approximated as a low-order polyno¬ 
mial, log(AF^™*) = y;. ai(log A)*, with a maximum de¬ 
gree ranging from two to four, and represents the stellar 
photospheric contribution by C Oph, as well as any un¬ 
known featureless dust extinction. The extinction optical 
depth is modeled using a linear combination of amor¬ 
phous silicates and H 2 O ice absorption spectra. We as¬ 
sume that these individual absorbing components have 
fixed profiles, but their relative contribution to the total 
optical depth may vary. 

We adopt amorphous silicate optical constants with 
stoichiometric compositions from th e melting and 
quenching experimental studies of iDorschner et al.1 
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Table 1 

Summary of Best-fit Optical Depth Parameters Toward Oph 


Poly-'‘ "fsil "fsil Tsil -fsil TH 2 O '^Fe 304 ’" T'pom'' xl 

Degree (total) (Mg 2 Si 04 ) (MgFeSi 04 ) (MgSiOa) 

Spitzer-IRS spectral decomposition models 

2 . 0.032 ± 0.009 0.014 ± 0.006 0.010 ± 0.006 0.009 ± 0.004 <0.006 <0.002 1.21 327 

3 . 0.031 ± 0.011 0.013 ± 0.007 0.011 ± 0.007 0.008 ± 0.004 <0.006 <0.002 1.21 326 

4 . 0.031 ± 0.011 0.013 ± 0.007 0.011 ± 0.007 0.008 ± 0.004 ■■■ <0.006 <0.002 1.21 325 

ISO-SWS spectral decomposition models 

2 . ■■■ ■■■ <0.008“^ ■■■ . 

2 . ■■■ ■■■ <0.009® ■■■ 2.34 1957 


Note. — Uncertainties are based on statistical errors in the Spitzer-IKS and ISO-SWS spectra only. 

^ Degree of the polynomial flux continuum. 

^ Optical depth upper limit estimate of the 16 /rm Fe 304 feature. 

^ Optical depth upper limit estimate of the 8.9 pm POM feature. 

^ Optical depth upper limit estimate of the 3.1 pm H 2 O ice feature for sub-micron-sized mantled grains. 

® Derived from the 95% upper confidence limit of the 3.53 pm H 2 O ice feature for thick mantles with central radius of a = 2.8 pm. 


and iHenning fc Stognienkol (jl99fi[ ). In particu¬ 
lar, we use Fe-free and Fe-bearing amorphous silicates 
with pyroxene (MgSiOs and MgFeSi206) and olivine 
(Mg 2 Si 04 and MgFeSi 0 4 ) compos it ions. Following the 
statistical approach of iMin et all 1)20031) . the optical 
properties of the amorphous silicate species are mod¬ 
eled using a distribution of hollow spheres (DHS) in 
the Rayleigh limit (A ^ 27ra, where a is the radius of 
the largest grains). In DHS theory, the optical proper¬ 
ties of hollow spherical particles are uniformly averaged 
from 0 < / < /max, where / is the fraction of the total 
volume occupied by the central vacuum inclusion. The 
value of /max indicates the degree of irregularity of the 
particles, with small values of /max representing nearly 
spherical particles. Adopting the maximum volume frac¬ 
tion o f /max = 0.7 d etermined for the Sgr A* line of 
sight (|Min et al.ll2007D . the resulting ensemble of parti¬ 
cles is representative of small, irregularly-shaped amor¬ 
phous silicate grains. 

For the H 2 O ice component, we adopt optical con¬ 
stants derived from the laborato ry spectrum of pure , 
amorphous H 2 O ice (T = 10 K: iHudgins et al.lll993D . 
Because irregularly shaped particle models have been 
largely successful at simul ating the profile shape o f in¬ 
terstellar ice features (e.g.- lPontoppidan et al.l[2?if)^ . the 
optical properties of H 2 O ice are modeled using a con- 
tinuo us distribution of ellipsoid s (CDE) in the Rayleigh 
limit (iBohren &: Huffmanlll983l) . However, we note that 
the calculated absorption cross section of H 2 O ice is not 
very sensitive to the choice of the adopted particle shape 
model (i.e., only small differences in the profile shapes 
of the absorption features are found if we instead adopt 
small, homogenous spherical mantles). 

Utilizing the model described in Equation © and 
the aforementioned solid-state absorption spectra, we 
employ a nonlin ear least-squares fitting technique 
(jMarkwardd 120091 ) to simulate the mid-infrared SED of 
C Oph. When present, blended silicate and H 2 O ice ab¬ 
sorption features occupy nearly all of the mid-infrared 
spectral region, and for this reason the Spitzer-IRS spec¬ 
trum is modeled over the entire 5.3-34 /rm wavelength 
range. In the fitting procedure, the initial values of the 
polynomial coefficients of the flux continuum and the 


scaling parameter of each absorbing component are ran¬ 
domized, with parameter convergence typically occurring 
after 30 or fewer iterations. Adopting a detection limit 
criterion of >1(T, the best-fit model parameters are eval¬ 
uated for significance upon convergence. If absorbing 
components are detected with a significance less than 
1(7, their corresponding scaling parameters are set to 
zero and the fitting procedure is repeated until all re¬ 
maining components satisfy the detection criterion. The 
formal uncertainties reported by the fitting procedure 
are derived from the statistical uncertainties within the 
Spitzer-IRS spectrum. If the resultant model is of good 
quality (i.e., the chi-squared per degree of freedom, xii 
is of order unity), the estimated parameter uncertainties 
are subsequently scaled by the square root of xt- 

The best-fit spectral decomposition model results to¬ 
ward C Oph are presented in Figure [2] for a second- and 
fourth-order polynomial flux continuum. The Spitzer- 
IRS spectrum and best-fit model spectrum are shown in 
the top panels, while the contributions by sightline ab¬ 
sorbers to the observed optical depth spectrum are pre¬ 
sented in the bottom panels. The best-fit optical depth 
parameters are listed in Table[T]for models using second-, 
third-, and fourth-order polynomial flux continua. 

3.1. Amorphous Silicate Composition 

In general, we find that the 5.3-34 pm optical depth 
spectrum toward C, Oph is adequately (y^ = 1.21) simu¬ 
lated by a linear combination of amorphous silicates with 
Mg 2 Si 04 , MgFeSi 04 , and MgSiOa compositions. Specif¬ 
ically, we find identical results for spectral decomposition 
models utilizing a third- or fourth-order polynomial flux 
continuum. In these cases, the Fe-free and Fe-bearing 
amorphous silicate species are detected with a signifi¬ 
cance ranging from 1.6tT-2(T. Furthermore, we find that a 
fit of equal quality may be obtained for a spectral decom¬ 
position model consisting of a second-order polynomial 
flux continuum. In this case, the Fe-free and Fe-bearing 
amorphous silicate components are detected at a simi¬ 
lar signihcance level ranging from 1.7a-2.3cr. The con¬ 
sistency between the best-fit optical depth results sug¬ 
gests that the amorphous silicate grain composition is 
not sensitive to the adopted polynomial flux continuum. 
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Finally, we find no evidence of Fe-bearing amorphous sil¬ 
icates with pyroxene composition (MgFeSi 206 ) or H 2 O 
ice along the line of sight. 

Depending on the choice of the adopted polynomial 
flux continuum, our analysis indicates that the peak sili¬ 
cate optical depth toward ( Oph is Tsu = 0.031 ± 0.011 or 
0.032 ± 0.009. In the local-diffuse ISM, the peak silicate 
optical depth is correlated with the ext inction along the 
line o f sight according to Ay /tsu ss 18.5 (|Roche fc AitkenI 
I1984H . Adopting a total-to-selective extinction ratio of 
Ry = Ay/Es-y = 2.55 ± 0.24 and a color excess of 
Es-y = 0.32 ± 0.04 from iValencic et id] (j2004H . we pre¬ 
dict a peak silicate optical depth of Tsii = 0.044 ± 0.007 
for the C Oph sightline, which is a factor of ^1.4 greater 
than our measured values. However, given the uncertain¬ 
ties associated with the measured and predicted values, 
we conclude that reasonable agreement is found among 
these quantities. If the small difference between the mea¬ 
sured and predicted values is real, then it might be ex¬ 
plained by partial masking of the observed silicate ab¬ 
sorption by emission from a warm dust component near 
COph. 
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3.1.1. Degeneracy Test for Silicate Components 

The robustness of the best-fit results was investigated 
by randomizing the initial values of the parameters 10^ 
times. On each occasion, we find that the spectral de¬ 
composition models converge to the same results. To 
further explore the possibility of degeneracy among the 
best-fit amorphous silicate components, the correlation 
coefficient, r, for each off-diagonal element is computed 
using 


where Ui and Uj are the formal parameter uncertainties 
of the ith and jth components, respectively, and covy 
is the covariance between the component pairs. Highly 
degenerate component pairs are anti-correlated and pos¬ 
sess a correlation coefficient near r = —1. Similarly, 
non-degenerate component pairs are uncorrelated and 
exhibit a correlation coefficient near r = 0. The most 
significant degeneracy is found between the Mg 2 Si 04 
and MgFeSi 04 components of the third- and fourth- 
order polynomial spectral decomposition models. For 
this component pair, we calculate a correlation coefficient 
of r = —0.84. Conversely, we find that only small degen¬ 
eracies, corresponding to r = —0.16 and r = —0.32, exist 
between the MgSi 03 :Mg 2 Si 04 and MgSi 03 :MgFeSi 04 
component pairs, respectively. 

In order to visualize the degeneracy between the Fe- 
free and Fe-bearing amorphous silicate components with 
olivine stoichometry, we recompute the fourth-order 
polynomial spectral decomposition model for 33 x 33 
combinations of the Mg 2 Si 04 and MgFeSi 04 compo¬ 
nents. All model parameters are held at their best-fit 
values except for the Mg 2 Si 04 and MgFeSi 04 parame¬ 
ters, which are varied over an optical depth range of 0 
< T < 0.032. Contours of the two parameter xt space 
are presented in Figure |3] and indicate the following ob¬ 
servables: (1) the best-fit {xt = 1-21) component pair 
from Table [T] is situated within the innermost contour, 
which is representative of the true global minimum. (2) 
The xi minimum is characterized by an ellipse with a 


Figure 3. Reduced chi-squared (x^) space of the fourth-order 
polynomial spectral decomposition model results. Contours of xi 
for the r(MgFeSi 04 ) and r(Mg 2 Si 04 ) parameters are illustrated 
for values ranging from 1.21 to 3.50. The shape of the Xu uiini- 
mum is characterized by an ellipse with a slope of approximately 
—0.8, indicating that the linear combination of the MgFeSi 04 and 
Mg 2 Si 04 components is well determined. Additionally, the abun¬ 
dance of elemental Fe (A^(Fe)/A’(H); black horizontal lines) con¬ 
tained in amorphous silicate grains is also illustrated for each model 
(see Section [4] for details of the calculation). Finally, the best-fit 
(x^ = 1.21) optical depth values and formal uncertainties from 
Tabled (white circle) are indicated for reference. 

slope of approximately —0.8. (3) The orientation of the 
xt minimum demonstrates that a linear combination 
of both components is required to accurately simulate 
the absorption profile by amorphous silicates toward C 
Oph. (4) Although the optical depth of the Mg 2 Si 04 and 
MgFeSi 04 components cannot be precisely disentangled 
from the total line of sight contribution, an indication of 
their dependence on one another is well determined. As 
such, increasing the optical depth of one component and 
decreasing the other yields a goodness-of-fit of similar 
quality {xt = 1-22). Hence, the two parameter xt space 
is instrumental in revealing both qualitative and quanti¬ 
tative information that could not otherwise be extracted 
from the formal Icr parameter uncertainties. 

3.2. H 2 O Ice Mantles 
3.2.1. Sub-Micron-Sized Grains 

Our model simulations of the Spitzer-IRS spectrum 
toward C Oph do not require a contribution from the 
12.6 ^m H 2 O ice libration feature. However, because this 
feature is blended with the strong vibrational modes of 
amorphous silicates, it may be difficult to detect along 
lines of sight where H 2 O ice is not abundant. In con¬ 
trast, the 3.1 ^m O—H stretching mode is intrinsically 
stronger than the libration mode (t 3 ,i/ti 2,6 ~ 4.5) and 
far less affected by overlap with other absorption fea¬ 
tures, and thus is a more reliable and sensitive tracer of 
H 2 O ice mantles on sub-micron-sized grains. 

Following a similar methodology as before, a spectral 
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Figure 4. Extinction cross sections for grain size distributions 
of thick H 2 O ice mantles on large grains. The weighted sum pro¬ 
files are presented for central grain radii of a = 2.0 (black line), 
2.2 (green line), 2.4 (blue line), 2.6 (red line), 2.8 (magenta line), 
3.0 (cyan line), and 3.2 pm (violet line). The individual scatter¬ 
ing (dotted violet line) and absorption (dashed violet line) cross 
sections are indicated for a grain size distribution with a central 
radius of a = 3.2 pm. 

decomposition analysis of the ISO-SWS spectrum was 
performed over the 2.8-3 .8 /im wavelength range using 
only the calculated absorption cross sections of H 2 O ice 
from Section [3l However, we find no evidence of the 
3.1 ^m H 2 O ice feature at the > 1(7 significance level. 
Consequently, from the standard error estimate between 
the observed optical depth spectrum and the second- 
order polynomial spectral decomposition model over the 
2.8-3.3 ^m spectral region, we deduce a limiting opti¬ 
cal depth of r < 0.008 for the 3.1 fim H 2 O ice fea¬ 
ture. We note that this result corresponds to a limit¬ 
ing optical depth of r < 0.002 for the 12.6 ^m H 2 O 
ice libration mode, which is a factor of ^2-3 less than 
the uncertainties reported for best-fit optical depth val¬ 
ues of the amorphous silicate components. Furthermore, 
the non-detection of H 2 O ice toward ( Oph is in ac¬ 
cord with the finding that surfaces of small (sub-micron- 
sized) grains begin to accumulate H 2 O ice mantles at 
column densities corresponding to silicate optical depths 
of Tsii 0.15 ± 0.03 towy d several prototypical dark 
clouds (|Whittet et al.l[20T^ and references therein). 

3.2.2. Very Large Grains 

Having found no convincing evidence of sub-micron¬ 
sized icy grains, we next consider the possibility of thick 
H 2 O ice mantles on very large grains (a > 1 /rm). The ex- 
ist ence of mic r on-siz ed interstellar grains was proposed 
bv iGreenberd ()1974[ 1. More recently, possible evidence 
for the presence of micron-s iz ed int erstellar grains have 
been found by iWang et al.l (|2014D . who modeled the 
“flat” mid-infrared extinction toward a variety of inter¬ 
stellar environments using a population of micron-sized 
grains. Furthermore, the recent analysis of three pos- 



Figure 5. ISO-SWS 2.4-3.9 /im spectrum and best-fit second- 
order polynomial spectral decomposition model for the C, Oph sight¬ 
line. Top panel: observed spectrum (black line), model fiux contin¬ 
uum (gray line), and best-fit model spectrum (blue line). Bottom 
panel: observed optical depth spectrum (black line) and best-fit 
model optical depth (blue line). The total extinction is modeled 
by a grain size distribution of thick H 2 O ice mantles with a central 
radius of a = 2.8 /im. 

sibly interstellar dust grains returned by the Stardust 
spacecraft provides tentative evi dence for the presence of 
micron-sized grains in the ISM ([Westphal et al.l[2014ll . 

Adopting amorphous silic a.te and H 9 O ice (T = 
25 K) optical c o nstan ts from iDraine fc Led (|1984f) and 
iMastrana et al.l (|2009H . respectively, the optical proper¬ 
ties of large, ice-coated spheric :al grains are modeled us- 
ing an extension of Mie theory (jBohren fc Huffmanlll983|l 
for particle radii in the range of 1 ^m < a < 4 /rm. We 
assume that one-eighth of the total particle volume is oc¬ 
cupied by the amorphous silicate core. In this thick man¬ 
tle limit, we note that the ice-coated extinction profile is 
similar to that of pure homogeneous H 2 O ice spheres. 
Finally, to smooth oscillatory structure in the calculated 
absorption cross sections, we sum over three grain radii 
of a — 0.2 /im, a, and a + 0.2 /im, with relative weights 
of 25%, 50%, and 25%, respectively. 

The extinction cross sections for grain size distribu¬ 
tions of thick H 2 O ice mantles are presented in Figure |4] 
with central radii of a = 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, and 
3.2 /im. The contributions from absorption and scatter¬ 
ing are separately shown for the largest grain size distri¬ 
bution, and demonstrate that the extinction cross sec- 
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Table 2 

Solid-State Column Densities Toward ^ Oph 


Poly.'^ S(sil) S(Mg 2 Si 04 ) S(MgFeSi 04 ) S(MgSi 03 ) S(Fe 304 )‘’ Af(H20) Ar(POM)‘= X(oli)d 

Degree (10“® g cm~^) (10“® g cm~^) (10~® g cm“^) (10“® g cm~^) (10~® g cm~^) (10^® cm“^) (10^® cm~^) (%) 


Spiizer-IRS spectral decomposition models 

2 . 1.1 ± 0.3 5.6 ± 2.4 3.8 ± 2.3 1.9 ± 0.8 <2.5 ■■■ <9.7 85 

3 . 1.1 ± 0.4 5.2 ± 2.8 4.2 ± 2.7 1.6 ± 0.8 <2.5 ■■■ <9.7 85 

4 . 1.1 ± 0.4 5.2 ± 2.8 4.2 ± 2.7 1.6 ± 0.8 <2.5 ■■■ <9.7 85 


ISO-SWS spectral decomposition models 

2 . ■■■ ■■■ ■■■ < 1 . 2 ® 

2 . ■■■ ■■■ ■■■ <3.2f 


Note. — Uncertainties are based on statistical errors in the Spitzer-IRS and 7SO-SWS spectra only. 

^ Degree of the polynomial flux continuum. 

® Based on the 16 pra Fe 304 optical depth. 

® Based on the 8.9 /rm POM optical depth. 

Mass fraction of amorphous silicates with olivine composition. 

® Based on the 3.1 /rm H 2 O ice optical depth of sub-micron-sized mantled grains. 

^ Based on the 95% upper confidence limit of the 3.53 /rm H 2 O ice optical depth for thick mantles with central radius of a = 2.8 /rm. 


tion of large grains are dominated by scattering at wave¬ 
lengths shortward of 2.8 71 m and longward of 3.3 71 m. 
Moreover, we find that the extinction profile shifts long- 
ward and decreases in amplitude with increasing grain 
radius. 

Employing a similar approach as before, a spectral 
decomposition analysis of the ISO-SWS spectrum was 
performed over the 2.5-3.9 /im wavelength range using 
only the aforementioned extinction spectra of very large, 
spherical H 2 O ice mantles. In general, we find that 
the ISO-SWS spectrum can be adequately {Xv = 2.34) 
simulated using a second-order polynomial flux contin¬ 
uum and grain size distributions with central radii in 
the range of 2 ^m < a < 3.2 ^m. The most signif¬ 
icant {3a) evidence of large H 2 O ice mantles is found 
for central radii of a = 2.6 and 2.8 71 m, which result in 
H 2 O ice optical depths of r = 0.006 ± 0.002 at 3.46 and 
3.53 fim, respectively. The best-fit spectral decomposi¬ 
tion model result is shown in Figure [5] for a grain size 
distribution with a central radius of a = 2.8 fim. We 
note that the 95% upper confidence limit of r < 0.009 
would only weakly affect the extinction law near 3.5 /rm 
(i.e., AAx/Eb-v = 1-086 x (0.009/0.32) = 0.03). This 
result demonstrates that the presence of H 2 O ice in the 
diffuse ISM may be revealed through scattering by thick 
icy mantles on large grains. The notion that such mate¬ 
rial might form within dense molecular clouds and sur¬ 
vive subsequent dispersal to more diffuse regions of the 
ISM may be supported by the longer lifetimes predicted 
for thick H 2 O ice mantles fe.g.. lOberg et al.ll2009[ l. 

4. SOLID-STATE ELEMENTAL ABUNDANCES 

The results from the spectral decomposition analysis 
allow constraints to be placed on the abundances of ele¬ 
ments contained within the relevant solids toward (/ Oph 
(i.e., those with vibrational modes in the spectral range 
covered by the available observations). A comparison be¬ 
tween gas-phase and solid-state abundances in this line 
of sight may be informative, particularly in regard to the 
question of any potential “missing” reservoir of interstel¬ 
lar O, as discussed in Section [T] 

In order to estimate the solid-state elemental abun¬ 


dances of Mg, Si, O, and Fe relative to H, we adopt a 
total hydrogen column density of N(R) = 1.42 x 10^^ 
cm~^ for the diffuse clouds toward ( Oph (iSavage et al.l 
[Hzi. However, the gas-phase interstellar absorption 
along this line of sight is known to origin ate from two 
separate diffuse cloud components (e.g., iHobbd IT^ 
iSnow fc MeversI I1979I ISavage et al.l 119921 1. The high- 
velocity component (vq = —27 km s“^) samples only the 
modest depletion associated with lower density clouds. 
In contrast, the low-velocity component (uq = —15 km 
s“^) samples the highly depleted, denser cloud medium, 
and account s for most of th e observed hydrogen along the 
line of sight (|,Ienkinsll2009l l. Adopting a synthetic hydro¬ 
gen column density of N(H) = 3.39 x 1 0 ^^ cm ~^ for the 
high-velocity component from iJenkinsI (|2009l l. we esti¬ 
mate a hydrogen column density of A^(H) = 1.39 x 10^^ 
cm“^ for the low-velocity component, which accounts 
for 97.6% of the total hydrogen column density toward 
COph. 

4.1. Amorphous Silicates 

Under the assumption that the absorption features in 
the Spitzer-IRS spectrum are entirely attributed to dust 
located in the highly depleted, low-velocity component, 
the relative abundance of element x contai ned in silicates 
is dete rmined following the prescription of lWhittet et al .1 

(IT^ : 

Njx) _ 1 fx,iTx^ , , 

7V(H) N{R) ^ 

where iV(H) = 1.39 x 10^^ cm“^, is the mass frac¬ 
tion of element x in the ith silicate component, fix is 
the atomic mass of element x, T\^i is the observed peak 
optical depth of the ith silicate component, and 

(«,.) = (4) 

^iPi 

is the shape-averaged peak mass absorption coeffi¬ 
cient for the ith silicate component with a shape- 
averaged absorption cross section per unit volume 














































Table 3 

Solid-State Elemental Abundances Toward f Oph 




Amorphous Silicates 


Fe304 

H 2 O Ice Mantles^ 

POM 

Poly.^ 

Degree 

JV{Mg) 

V(H) 

(ppm) 

N(Si) N(0) 

]V(H) iV(H) 

(ppm) (ppm) 

iV(Fe) 

vW 

(ppm) 

N{0) iV(Fe) 

JV(H) JV(H) 

(ppm) (ppm) 

N{0) N{0) 

Ar(H) JV(H) 

(ppm) (ppm) 

JV(0) 

vTiiy 

(ppm) 


Spiteer-IRS spectral decomposition models 


2 . 

. 52 ± 16 

34 ± 10 

130 ± 39 

10 ± 6 

<19 

<14 

<7 

3 . 

. 49 ± 19 

33 ± 11 

126 ± 45 

10 ± 7 

<19 

<14 

<7 

4 . 

. 49 ± 19 

33 ± 11 

126 ± 45 

10 ± 7 

<19 

<14 

<7 


/50-SWS spectral decomposition models 


2 


<9 <23 


Elemental abundances absent from the gas-phase'^ 


Protosolar^ 

42 ± 4 

34 ± 3 

230 ± 11 

35 ± 3 

Protosolar-l-GCE® 

46 ± 1 

41 ± 4 

282 ± 

48 ± 1 

F and G stars^ 

42 ± 

39 ± 

282 ± 207 

41 ±32 

B stars® 

34 ± 4 

30 ± 4 

268 ± 72 

33 ± 2 


Note. — Uncertainties are based on statistical errors in the Spitzer-IRS and ISO-SWS spectra only, 
assuming a hydrogen column density of Af(H) = 1.39 xlO^^ cm~^ for the low-velocity {vq = —15 km s“^) 


Abundances are calculated 
component. 


^ Degree of the polynomial flux continuum. 

^ Abundance of elemental O in sub-micron-sized and very large H 2 O ice mantles, respectively. The latter is derived from the 95% upper 

confidence limit of the 3.53 fim H 2 O ice optical depth for thick mantles with central radius of a = 2.8 fim. _ _ 

^ Based on the difference between reference abundances and the observed gas-phase abundances derived from [Jenkins! II2009I) : 1.9 ± 
0.1 ppm (Mg), 1.6 ± 0.1 ppm (Si), 307 ± 30 ppm (O), and 0.18 ± 0.01 ppm (Fe). 

JAs£luiM_e^^ll 1200^ 

® [XsphmT'er'arf 11255^ ; IChiappini et al] II2003I 1 
^ IBensbv et al.l 11200511 1 Tbodders et al.l 1200911 
^ INieva &: Przvbillal 120121 1 


{Cabs,i)/Vi and specific density pi. Using the laboratory- 
measured specific densities oi p = 3.2, 3.71, and 
2.71 K cm~^ for Mg 2 Si 04 , MKFeSi 04 , and Mg SiO.'^ 
(jHcnning fc StognienkolfTo^ iDorschner et idlllQQhl) . we 

calculate shape-averaged peak mass absorption coeffi¬ 
cients of (ka) = 2510, 2624, and 4853 cm^ g“^, respec¬ 
tively. 

The total and individual mass column densities, E = 
T\/{n\), of amorphous silicates, as well as the mass 
fraction of amorphous silicates with olivine composi¬ 
tion, are presented in Table [2] The solid-state ele¬ 
mental abundances of Mg, Si, O, and Fe are listed 
in Table |3l For comparison, reference elemental dust 
abundances, based on average abundances in young F 
and G stars (jBensbv et al.l 1200 51: iLodd ers et al.l I2OO90 . 
protosolar abundances (jAsnlund et al.l l2009| ). proto¬ 
solar abundances with Galactic chemical enrichment 
(GGE) over the past 4.56 Gyr (lAsolund et al.l 120091 : 
iChiaDpini et al.l 1200311. and a verage abundances in B 
stars (|Nieva fc Przvbillal[2012t) . are calculated using pre¬ 
viously published and revised gas-phase column densi- 
ties of the low -velocity cloud component toward C, Oph 
(|Jenkinsll2009L and references therein). 

In general, we find that solid-state elemental abun¬ 
dances derived from the spectral decomposition models 
are in good agreement with those predicted from stellar 
and solar abundances. Specifically, we find that the ele¬ 
mental abundance of Mg is slightly-to-moderately over¬ 
consumed, requiring ^7%-53% more than the available 
abundance. Similarly, our model results indicate that 
nearly all (^80%-100%) of the available elemental Si 
abundance is present in amorphous silicates. (When con¬ 


sidering the average reference abundance in B stars, we 
note that the elemental abundance of Si is slightly over¬ 
consumed by ^10%-13%.) In contrast to the substantial 
fraction of elemental Mg and Si contained in amorphous 
silicates, we find that ^70%-79% of the available ele¬ 
mental Fe abundance is expected to be present in some 
other solid-state species. Finally, our model results in¬ 
dicate that only ^45%-57% of the elemental O that is 
missing from the gas-phase is present in amorphous sili¬ 
cate grains. Independent of the adopted reference abun¬ 
dances, these results suggest that additional reservoirs of 
interstellar Fe and O are needed along the line of sight 
toward C Oph. 

4.2. Additional Repositories of Interstellar Fe 
4.2.1. Fe-rich Amorphous Silicates 

Similar to prev i ous diffuse ISM studies (e.g., 
Keinoer et al.l 12004 iChiar fc TielensI 120061 I Min et al] 
200711 . Fe-bearing amorphous silicates with equal 
amounts of elemental Mg and Fe have been utilized to 
simulate the absorption profile of amorphous silicates to¬ 
ward f Oph. To investigate whether a significant frac¬ 
tion of interstellar Fe may be present in amorphous sil¬ 
icates with Fe-to-Mg ratios greater than unity (Fe/Mg 
> 1), the Spitzer-IRS spectrum toward ( Oph is further 
modeled using the Mg 4 Fe 6 Si 5 02 o species instead of the 
MgFeSi 04 component. However, we find that the slightly 
Fe-rich component is not utilized by the nonlinear least- 
squares fitting routine; this is true for all (2nd-, 3rd-, and 
4th-order) polynomial flux continua considered in Sec¬ 
tion |3l This result is likely a consequence of the broad¬ 
ened spectral profiles possessed by Fe-rich amorphous sil- 
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Figure 6. Spitzer-IRS spectrum and best-fit fourth-order poly¬ 
nomial spectral decomposition model for the Oph sightline. Top 
panel: observed spectrum (black line), model flux continuum (gray 
line), and best-fit model spectrum with (solid blue line) and with¬ 
out (dotted blue line) Fe 304 . Bottom panel: observed optical 
depth spectrum (black line) and best-fit model optical depth with 
(solid blue line) and without (dotted blue line) Fe 304 . The to¬ 
tal absorption is modeled as a linear combination of Fe 304 (ma¬ 
genta line) and amorphous silicates with Mg 2 Si 04 (green line), 
MgFeSi 04 (violet line), and MgSi 03 (red line) stoichiometries. A 
limiting optical depth of r < 0.006 for the 16 ^m Fe 3 04 feature is 
assumed. 

icates (|Dorschner et al.iri995D . and further demonstrates 
that stoichiometries with Fe/Mg < 1 are required when 
simulating the absorption profile of diffuse ISM silicates. 
Finally, if we assume the unlikely notion that MgFeSi 04 
is the only olivine-like species present toward C Oph, we 
find that the abundance of elemental Fe is increased by 
only ^10 ppm, as indicated in Figure [31 Therefore, we 
conclude that some other Fe-bearing compound must ac¬ 
count for the ^^23-38 ppm of elemental Fe missing along 
the line of sight. 

4.2.2. Fe Oxides and Metallic Fe 

In addition to its presence in silicate grains, in¬ 
terstellar Fe could exist in pure metallic Fe or Fe- 
beari ng oxide particles of Fe O, 7 -Fe 203 , and Fe 304 
(e.g., iDraine &: Hensle"^ I2013L and references therein). 
Metallic Fe possesses no active vibrational modes in the 
mid-infrared wavelength region and contributes only to 
the continuum opacity of dust. Pure metallic Fe par¬ 
ticles have been found toward the circumstellar dust 


shells surrounding a symptotic giant branch stars (e.g., 
iKemoer et al.l 1200211 . but may be unlikely to remain 
pristi ne in the ox ygen-rich environment of the diffuse 
ISM (iJonesI irpOOH . In contrast, Fe-bearing oxides are 
considered potent ial constituents of interstellar grains 
(ISofia et al.l II994I1 . Depending on the precise compo¬ 
sition of these species, Fe-bearing oxides exhibit single 
or multiple vibr ational resonances within the 17-28 /nn 
spectral region (iHenning et al.lll995l : iGlotch et al.l 1200^ 

IGlotch fc RossmanI 1200911 . However, because their spec¬ 
tral features are blended with the strong bending mode 
of amorphous silicates, Fe-bearing oxides have not been 
detected in the diffuse ISM, and thus are not usually con- 
sidere d to be a significant co mponent of interste llar dust 
(e.g., IChiar fc TidensI 1200611 . More recently, iLi et al.l 
(1201311 found that nanoparticles of FeO are not likely 
the carrier species responsible for the 21 /rm emission 
observed toward post-asymptotic giant branch stars. 

To examine whether Fe-bearing oxides are a major 
reservoir of elemental Fe toward ^ Oph, we consider ir¬ 
regularly shaped particles of FeO (wiistite) and Fe 304 
(magnetite) as potential line of sight absorbers. While 
7 -Fe 203 (maghemite) is known to possesses broad reso¬ 
nance features near 17 and 27 /iin, its absorption opacity 
contributes primarily to the featureless dust continuum, 
and a non-detection of the features does not provide a 
useful upper limit estimate on the abundance of 7 -Fe 203 . 
For this reason, we do not include 7 -Fe 203 in our spectral 
analysis. 

As a result of the low signal-to-noise ratio within the 
19-35 ^m region of the 5piteer-IRS optical depth spec¬ 
trum, we are unable to accurately quantify the pres¬ 
ence of Fe-bearing oxides along the sightline using the 
spectral decomposition analysis from Section [3] Alter¬ 
natively, from the standard error estimate between the 
observed and modeled optical depth spectrum within the 
15-30 /im wavelength region, we deduce a limiting optical 
depth of T < 0.009 for the 22 /im FeO fe ature. Assum¬ 
ing a specific density oi p = 5.70 g cm“^ ([Henning et al.l 
119951) and a DHS (/max = 0.7) shape-averaged peak ab¬ 
sorption cross section of {Cahs)/V = 13590 cm“^, we set 
a limit of <23 ppm for the abundance of elemental Fe 
in FeO grains. (In this context, a limiting abundance of 
<23 ppm for elemental O in FeO is also deduced.) Nev¬ 
ertheless, we emphasize that the generous limit on the 
optical depth of FeO, which is likely a consequence of 
the appreciable noise present in the 19-30 pm spectral 
region, would considerably affect the shape of the optical 
depth spectrum toward / Oph, and therefore should be 
considered with caution. 

Similarly, from the standard error estimate between 
the observed and modeled optical depth spectrum within 
the 13-18 pm wavelength region, we calculate a limit¬ 
ing optical depth of r < 0.006 for the 16 pm Fe 304 
feature. In contrast to the FeO investigation, we note 
that the limit on the optical depth of Fe 3 04 does not 
strongly affect the appearance of the 13-30 pm re¬ 
gion of the best-fit models from Section 13.11 as illus¬ 
trated in Fig ure El Adopting a spec ific density of p = 
5.18 g cm“^ (IDraine &: Hensle)^l2013ll and a DHS shape- 
averaged peak absorption cross section of {Csbs)/V = 
12260 cm“^, we deduce limiting abundances of <14 ppm 
and <19 ppm for elemental Fe and O contained in Fe 304 , 
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respectively. In summary, these results establish that 
Fe 304 may account for <37-61% of the missing elemen¬ 
tal Fe abundance toward ( Oph. 

4.3. Missing Reservoirs of Interstellar 0 

As previously noted, the rate at which the observed 
depletion of elemental O varies with mean line of sight 
density in the diffuse ISM is more rapid than can be 
accounted for by exchange of atoms with amorphou s sil¬ 
icates and metallic oxide particles (jJenkinsI I2009D . In 
consequence, a substantial fraction of interstellar O is 
thought to reside in an additional reservoir, presumably 
in combination with some other abundant element (s). 
Given that hydrogen and carbon are the only reactive 
elements with sufficiently large interstellar abundance 
(and given the improbability of substantial amounts of 
O 2 existing in such environments), it follows that H 2 O 
ice and organic refractory compounds are attrac tive can¬ 
didates as carriers of the missing elemental O (|Jenkinsl 
[20091: IWhitt^[20lnll . 

4.3.1. H 2 O Ice 

The abundance of elemental O contained in small H 2 O 
ice mantles may be directly inferred from its line of sight 
column density (jAllamandola fc Sandfordl[l988D : 

iV(H20) (5) 

where A is the laboratory-derived intrinsic band strength 
and f Tc,dv is the integrated optical depth in the 
wavenumber domain, v. Adopting an i ntrinsic band 
stren gth of A = 2.0 x 10“^® cm molecule”^ ([Hagen et al.l 
[MID for the 3.1 /rm stretching mode, we calculate a 
limiting column density of iV(H 2 0) < 1.2 x 10^® cm“^. 
The abundance of elemental O contained in small H 2 O 
ice mantles toward the low-velocity component of C Oph 
is thus no more than 9 ppm (i.e., no more than 6%- 
9% of the missing elemental O abundance). We note 
that an even stricter limit (<2 ppm) was reported to¬ 
ward the hi ghly reddened (Ay R i 10) luminous star Cyg 
OB2 no. 12 ([Whittet et al.lll9^ . a sightline also thought 
to be dominated by di ffuse ISM material at compara¬ 
tively high density fe.g.. lGeballe et aOl999l : lGredel et al.l 
I200lr) . These results demonstrate that H 2 O ice mantles 
on sub-micron-sized grains are not a significant reservoir 
of elemental O near the diffuse-dense ISM transition. 

By analogy with Equation^ the abundance of elemen¬ 
tal O in thick H 2 O ice mantles on large grains is given 

by 

N{0) _ 1 Tx pVn^o /gN 

A(H) “ A^(H) (TgQext.A ISmn ’ ^ ’ 

where Cg is the geometric cross section of the grains 
(tto^), HH 2 O is the volume of the mantle, and Qext.A is 
the peak extinction efficiency of H 2 O ice. Assuming a 
spe cific density of p = 1.1 g cm“^ for amorphous H 2 O 
ice (INarten et al.l[l97^ . the abundance of elemental O 
in thick H 2 O ice mantles is presented in Figure [7] for 
grain size distributions with central radii in the range of 
2 /im < a < 3.2 /rm. Gonsidering the 95% upper confi¬ 
dence limit for a grain size distribution central radius of 
a = 2.8 /im, we find that as much as ~23 ppm of elemen¬ 
tal O may be stored in thick H 2 O ice mantles on large 



Figure 7. Abundance of elemental O that may be contained in 
thick H2O ice mantles on large grains along the ^ Oph sightline for 
grain size distributions with central radii of a = 2.0, 2.2, 2.4, 2.6, 
2.8, 3.0, and 3.2 pm. The magnitude of each data point represents 
the significance of the best-fit H2O ice optical depth parameters, 
which varies over the Icr-Sir level. 

grains toward f Oph. Although this result is more than 
a factor of two larger than the upper limit abundance 
obtained for small icy grains, it indicates that thick H 2 O 
mantles on large (a = 2.8 /tm) grains can only account 
for 15%-23% of the missing elemental O abundance. 

4.3.2. Organic Compounds 

We next consider organic refractory matter as a possi¬ 
ble repository of interstellar O toward C Oph. It is well- 
known that laboratory experiments simulating the ener¬ 
getic processing (e.g., ultraviolet photolysis or ion bom¬ 
bardment) of interstellar ices lead to the production of 
organ ic refractory matter with significant oxygen content 


fe.e:.. IGreenbere et al.l 119951: 

Pendleton & Allamandolal 

120021: IMuhoz Garo & Schuttel 

20031) . These oroducts in- 


elude species with carbonyl molecular groups (esters, 
amides, and carboxylic acid salts) and species related 
to polyoxymethylene (POM), a polycrystalline material 
consisting of linear chains of formaldehyde molecules 
[(—GH 2 O—)„]. The ester and amide species possess 
minor G=0 stretching modes near 5.7 and 6.0 /^m, 
respectively, while carboxylic acid salts are known to 
exhibit a strong GOO~ str etching mode at 6.3 fim 
([Munoz Garo fc Schuttdl2003fl . However, apart from the 
detection of weak ^^6 iim absorption toward carbon-rich 
Wolf- Rayet stars ([Schutte et ai.l 119981: iGhiar &: TielensI 
I2OOIII . evidence for the presence of these 0-bearing or¬ 
ganic compounds in the diffuse ISM has remained elusive. 

The quality of the 5-8 /im region of the Spitzer-IRS op¬ 
tical depth spectrum toward C Oph does not permit any 
meaningful profile comparison with laboratory spectra of 
carbonyl molecular groups. In consequence, we make no 
attempt to constrain the possible fraction of interstellar 
O contained in carbonyl groups, and further consider the 
possibility of POM-like species instead. 
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Table 4 

Summary of the Observed Elemental Budget Toward ^ Oph 


Material 

JV(Mg) 

TvThT 

JV(Si) 

jvTHy 

JV(0) 

nThT 

A^(Fe) 

IvW 

Reference 


(ppm) 

(ppm) 

(ppm) 

(ppm) 


Adopted total . 

47.9 ± 

42.7 ± 

589 ± 11 

47.9 ± 

1 

Gas . 

1.9 ± 0.1 

1.6 ± 0.1 

307 ± 30 

0.18 ± 0.01 

2 

Silicates. 

49 ± 19 

33 ± 11 

126 ± 45 

10 ± 7 

3 

Unidentified . 

-3 -h 20 

8 ± 12 

156 ± |0 

38 ± 8 


H 2 O Ice (sub-micron-sized) 



<9 


3 

H 2 O Ice (a 3 pm) . 



<23 


3 

Fe304 . 



<19 

<14 

3 

POM . 



<7 


3 

Other . 

-3 -h 20 
^ 19 

8 ± 12 

(98-156) ± |0 

(24-38) ± 8 



References. — fl BAsplund et ^ 1120091b FChiappini et all II2003I) : (2) IJenkinsI ll2009l'l : (3) this work. 
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Figure 8. Total elemental abundances along the Oph sightline 
with respect to reference abundances. The refere nce value s are 
based on average abundan ces in young F and G stars dBens^ 

2Q0f: flTodders et al.ll2003) . protosolar abundances dAsplund 
200£ ), protosolar abundances with Galactic chemical enrichment 

’200c), and average abun- 
201^). The gas-phase el- 


et al.l 

plund et al.l 


JAsplund et al.l 


dances in B stars ^^Jiev^^^PrzvEiHa 


Chiar )Dini et al. 


emental abundances of Mg, Si, O, and Fe (magenta) are based 
on measured column densities of the low-velocity (uq = —15 km 
s“^) component Iljenkindl2009l . and references therein). The solid- 
state abundances (blue) are representative of the total number of 
atoms contained in amorphous silicate grains, and are derived from 
the fourth-order polynomial spectral decomposition model. Upper 
limit abundance estimates are indicated for elements contained in 
H 2 O ice mantles (green), Fe 304 (violet), and POM (red). The 
former represents the abundance of elemental O that may be con¬ 
tained in H 2 O ice mantles on both sub-micron-sized (<9 ppm) and 
very large (<23 ppm) grains. 


To explore the presence of POM-like species in the dif¬ 
fuse ISM, we adopt optical consta nts derived from th e 
laboratory spectrum of pure POM (| Whiffet et al.lHOT^ . 
The optical properties of POM are simulated using the 
DHS grain shape model (/max = 0.7) in the Rayleigh 
limit. We find that the POM species possesses a strong 


C—O—C vibrational mode at 8.95 /im, as well as several 
weaker features at 6.82, 7.34, 8.05, 11.36, and 11.80 ^m. 
Following a similar methodology as before, a spectral de¬ 
composition analysis of the Spitzer-IRS spectrum is per¬ 
formed using the absorption spectra of amorphous sil¬ 
icates and POM. However, we find no evidence of the 
8.9 ^m POM feature at the >lcr significance level. Al¬ 
ternatively, from the standard error estimate between the 
observed and modeled optical depth spectrum over the 
8.5-9.1 /rm wavelength region, we calculate a limiting 
optical depth of r < 0.002 for the 8.9 /im POM feature. 
Adopting an in trinsic band strengt h of A = 9.7 x 10“^® 
cm molecule"^ (|Schutte et al.lll99^ . the limiting column 
density of POM toward this diffuse sightline is 7V(POM) 
< 9.7 X 10^® cm“^. Thus, we conclude that no more 
than <7 ppm of elemental O is expected to be found in 
POM along the C Oph line of sight. 

5. DISCUSSION 

We have shown that the 5-36 /im Spitzer-IRS spec¬ 
trum toward / Oph can be adequately simulated using an 
ensemble of small, irregularly shaped amorphous silicate 
grains of olivine (Mg 2 Si 04 and MgFeSi 04 ) and pyrox¬ 
ene (MgSiOa) compositions. The average stoichiometry 
of the grains corresponds to Mg/Si = 1.5, 0/Si = 3.8, 
and Fe/Si = 0.3, indicating that the amorphous silicate 
mass fraction toward the diffuse sightline is dominated by 
olivine-like compositions (Aqu ~ 85 %). We no t e that this 
result is in agreeme n t with those of I Min et al.l (j2007t) and 
Ivan Breemen et ^ (j201ir i. who found that amorphous 
silicates in the diffuse ISM are Mg-rich and primarily 
olivine-like in composition. 

Presolar silicate grains found in primitive meteorites 
possess both stoichiometric (i.e., pyroxene- and olivine¬ 
like stoichiometries) and non-stoichiome tric composi¬ 
tions , with a wide ran ge of Mg/Fe ratios (jNeruven et al.l 
l2010l : lBose et al.ir2012D . These grains, identified by virtue 
of isotopic anomalies, are “stardust”, condensed in stel¬ 
lar outflows. Approximately 20%-25 % of stardust grains 
found in meteorites are crystalline (jNguven et al.l 12(11(1 
and references therein), and are thought to have under¬ 
gone only minimal p rocessing during pa ssage through 
the diffuse ISM (e.g.. IWllmer et ^120091) . In contrast, 
only <2% of interstellar silicates are of crystalline na¬ 
ture (|Kemper et al.l[2005l) . implying that stardust grains 
comprise at most ^10% of interstellar silicate material 
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(lDrainell200^ . Nevertheless, the varied silicate compo¬ 
sitions found among stardust grains are not in conflict 
with our findings on the bulk composition of diffuse ISM 
silicates toward ( Oph. 

Upon further investigation, we find no conclusive evi¬ 
dence for the presence of Fe-bearing oxides, POM species, 
or H 2 O ice mantles on small grains along the line of 
sight. However, from the spectral decomposition analy¬ 
sis of the ISO-SWS spectrum, we find tentative evidence 
for the presence of thick H 2 O ice mantles on large grains 
(a = 2.8 fim). These results establish that amorphous 
silicate grains account for at least 126 ± 45 ppm of ele¬ 
mental O toward ( Oph, while Fe 304 , POM species, and 
H 2 O ice mantles on sub-micron-sized grains can only ac¬ 
count for <19 ppm, <7 ppm, and <9 ppm of elemental 
O, respectively, and therefore are not significant reposi¬ 
tories of interstellar O along this line of sight. Further¬ 
more, tentative evidence for thick H 2 O ice mantles on 
large grains suggests that <23 ppm of elemental O may 
be present in icy grain mantles with radii a = 2.8 ^m, 
thereby establishing that H 2 O ice mantles may collec¬ 
tively contain <32 ppm of elemental O. 

The total elemental abundance of Mg, Si, O, and Fe in 
the low-velocity component (uq = —15 km s“^) toward 
C Oph are presented in Figure [51 The gas-phase abun¬ 
dances are inferr ed from t he rev ised column density mea¬ 
surements from iJenkinsI (|2009l and references therein), 
while the solid-state elemental abundances are derived 
from the fourth-order polynomial spectral decomposi¬ 
tion model. For co mparison, average abundances in 
youn g F and G stars (iBensbv et al.ll200^ iLodders et al.l 
1200911 . protosolar abundances ( Asplund et al.ll2009ll pro- 
tosolar abundances w ith GCE ( Asplund et al.l l2009t 
iGhiappi ni et al.1 1200311. and average abundances in B 
stars ([Nieva fc Przvbillal 120121 1 are also provided. The 
elemental budget toward C Oph is summarized in Ta¬ 
ble 01 wh ere we adopt protosolar abundances augmen ted 
by GCE (| Asplund et al.l[200^ IGhiappini et ^1200311 as 
the best estimate for interstellar abundances in the solar 
neighborhood. 

For elemental Mg and Si, we calculate total abun¬ 
dances of fV(Mg)/7V(H) = 51 ± 19 ppm and fV(Si) /N (H) 
= 35 ± 11 ppm, respectively. These results are consistent 
with the adopted reference abundances, indicating that 
essentially all elemental Mg and Si along the line of sight 
are present in amorphous silicates. In contrast, the de¬ 
tected abundance of iV(Fe)/7V(H) = 10 ± 7 ppm suggests 
that ~38 ± 8 ppm of elemental Fe resides in compounds 
other than amorphous silicates. From our upper limit 
abundance estimate, we find that Fe 304 could account 
for <37% of the missing elemental Fe abundance. 

In the case of elemental O, we find that the observed to¬ 
tal (silicates plus gas) abundance of A^(0)/iV(H) = 433 ± 
54 ppm is largely inconsistent with the adopted reference 
abundances presented in FigurejH Under the assumption 
that the protosolar elemental O abundance with GGE 
{N{0)/N{H) = 589 ppm) is representative of the cur¬ 
rent, local diffuse ISM, then as much as ^^156 ppm of ele¬ 
mental O is unaccounted for along the line of sight toward 
C Oph. (Note that the same result is obtained if we in¬ 
stead adopt the average elemental O abundance in young 
F and G stars as the true interstellar standard.) Ou r con¬ 
clusions are in accord with those of iJenkini (l2009f) . who 
showed that spatial variations in the consumption of el¬ 


emental O over many lines of sight cannot be accounted 
for by exchange with amorphous silicates alone. More¬ 
over, our limits on the abundances of H 2 O ice, POM, 
and Fe 304 imply that these materials can account for 
at most ~37% of the missing ~156 ppm of elemental O, 
leaving ~98 ppm of elemental O unaccounted for along 
the line of sight. Thus, we have demonstrated that, other 
than absorption by amorphous silicates, no conclusive 
evidence for an abundant 0 -bearing solid exists in the 
Spitzer-IRS or ISO-SWS spectra toward C Oph. There¬ 
fore, if a missing compound of elemental O is present in 
the diffuse ISM, it must not contribute significantly to 
the extinction opacity in the near- to mid-infrared wave¬ 
length region. 

As discussed by IJenkinsI (|2009f) , hidden reservoirs of 
elemental O may reside on very large grains (a 3> 1 pm) 
that are nearly opaque to infrared radiation. However, 
as demonstrated in this study, the presence of thick 
H 2 O ice mantles on grains with radii as large as 3.2 pm 
would introduce conspicuous structure within the 2.5- 
4.5 ^m spectral region (see Figure 01 ). Moreover, we 
note that the observed interstellar extinction and po¬ 
larization from ultraviolet to near-infrared wavelengths 
toward ( Oph suggest no evidence for the presence of 
unusually large grains along the line of sight. The star 
possesses an unexceptional ultraviolet extinction curve, 
with a somewhat weaker-than-ave rage 2175 A extinction 
bump ([Fitzpatrick fc Massalll99^ . Furthermore, the ra¬ 
tio of total-to-selective extinction (Ry = Av/Eb-v ~ 
2.6) appears to be somewhat less than the diffuse-ISM 
mean fSection 13.11) . suggesting a preponderance of rel¬ 
atively small grains; although the wavelength of maxi- 
mum linear polarizatio n toward C Oph (Amax ~ 0.59 ^m; 
iSerkowski et al.1 1197511 is marginally greater than the 
diffuse-ISM median. Thus, if thick H 2 O mantles on large 
grains are present along the line of sight, they would con¬ 
tribute relatively little extinction at optical wavelengths, 
and very little reddening, and would be difficult to de¬ 
tect. 

The abundance of interstellar O in all constituents of 
the diffuse ISM (both solid and gaseous) may be ex¬ 
amined through the detailed structure of the 540 eV 
oxygen K-shell absor ption-edge fea t ure t o ward bright 


X-ray 

sources fe.g.. IPaerels et al.l 120011 iTakei et al. 

2002 

:ICostantini et alJl2n05HBaume:artner Kr. Mushotzkv 

2006 

). However, the largest grains (a > 1 ^m) become 


optically thick to X-rays at energies near 540 eV, result¬ 
ing in a nearly constant extinction cross section over the 
absorption-edge. Detailed spectral decomposition anal¬ 
ysis of this feature toward nine Galactic sightlines sug¬ 
gests that ~15%-25% of the total elemental O abundance 
may be present in solids such as amo rphous silicates and 
H 2 O ice mantles ([Pinto et al.ll201^ . Nevertheless, de¬ 
generacies between the absorption features produced by 
different species prevent the precise identification of the 
chemical composition of the 0-bearing reservoirs. Thus, 
spectral modeling of the oxygen K-shell absorption-edge 
offers another possibility of identifying significant repos¬ 
itories of elemental O on sub-micron-sized grains, but 
is unable to identify interstellar O that may be present 
in large grains (a > 1 ^m). In conjunction with detailed 
spectral decomposition analyses of currently available in¬ 
frared spectra, future spectroscopic observations with the 
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James Webb Space Telescope may provide the best op¬ 
portunity for detecting thick H 2 O ice mantles on large 
grains, if they exist. 
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